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SUMMARY 

Supersonic profiles of minimum pressure drag for a given thickness 
ratio and for a given area have been determined with the use of a non- 
linear pressure relation and are compared with minimum- drag profiles 
found by linearized theory. The results show that the profiles are 
determined with sufficient accuracy by linear theory over the entire 
supersonic Mach number range since the drag coefficients for these 
profiles are only slightly higher than those for optimum profiles deter- 
mined by no n li n ear theopy. Linear theory appears to be adequate for 
determining profiles of minimum drag for other auxiliary structural 
conditions since moderate deviations from the optimum shape have only 
a small influence on the pressure drag. 

The parameters determining the airfoil shape for a given thickness 
ratio found by both the linear and nonlinear theory are presented in 
graphs as a function of the base pressure coefficient. With the use of 
these results, the optimum profiles for any stream Mach number and thick- 
ness coefficient are readily determined. A comparison of the pressure 
drag coefficients for optimum profiles determined by linear and no n li n ear 
theory is presented for the Mach number range from 1.5 to 10.0. In 
addition, several optimum profiles for a given area have been calculated 
by both the linear and nonlinear theory. 


INTRODUCTION 


Drougge (reference l) has determined the airfoil section shape for 
minimum pressure drag at supersonic speeds subject to such auxiliary 
conditions as given bending and torsional stiffness. These calculations 
were made by using the linearized expression for the pressure coefficient^ 
the effect of a base was not considered. Recently, Chapman (reference 2) 
has shown that the section shape for minimum pressure drag as determined 
by linearized theory may have a blunt trailing edge. , The use of linear- 
ized theory for determining optimum profiles facilitates the mathematical 
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development; however, the results are subject to question particularly 
at high- supersonic Mach. numbers. 

The purpose of the present paper is to compare the section shapes 
for minimum pressure drag (subject to certain auxiliary conditions) 
determined by linear and nonlinear theory in order to estimate the 
errors introduced by the linearized form of the pressure coefficient 
and to determine its range of validity for calculations of this nature. 
For this purpose, it was considered sufficient to examine two problems. 
The problems chosen were the determination of the profile for minimum 
drag for a given thickness ratio and the determination of "the profile 
for minimum drag for a given area. 

The nonlinear form of the expression for the pressure coefficient 
used in the present analysis is derived in reference 3 where ix is 
shown to be in excellent agreement with the exact expression for stream 
Mach numbers greater than 1.5* The variation of base pressure coeffi- 
cient with stream Mach number was assumed in order to determine actual 
airfoil shapes. This base-pressure curve was based in part upon experi- 
mental data and the known variation of vacuum pressure coefficient with 
Mach number. 


NOTATION 


X chordwise distance from leading edge 

Y airfoil ordinate 

if 

c airfoil chord; also, 

7 ratio of specific heat at constant pressure to specific heat 

at constant volume 



M^, stream Mach number 

m = \JlA^ - 1 

P pressure coefficient 
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t thickness coefficient 



■5 flow deflection angle 

0 •= md 

c^ drag coefficient 

A airfoil area 

Subscripts: 

1 front surface 

2 rear surface 

b base 

v vacuum 

w refers to double -wedge profile 

opt optimum 

TE trailing edge 


ANALYSIS AND DISCUSSION 


The pressure drag of an airfoil is the drag due to the normal forces 
acting over the airfoil surface. For a profile with a blunt trailing 
edge the base drag must be added to the drag due to the normal forces 
over the forward surfaces to give the total pressure drag. With the 
surface pressure coefficient denoted by P, the airfoil shape by Y(X), 
the chord by c, the base height by and the base pressure coef- 

ficient by P^ the pressure drag coefficient c^ of a thin symmetric 
two-dimensional airfoil at zero angle of attack is given approximately by 



( p - p b) 
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X Y Y-u 

or with, x = — , y = — , and y^ = — , the pressure drag may 


UC W-L U.U OCU 



p S 331 - Vi. 


( 1 ) 


Profile for minimum drag for a given thickness ratio .- In order to 
determine the profile for minimum drag for a given thickness ratio, it 
is reguired to determine the function y(x) and the base height y-^ 
which make c^ a minimum. For supersonic flow, to a high order of 

approximation, the pressure coefficient is a function only of the local 
dy 

slope and the stream Mach number. . By use of the calculus of varia- 


tions, any integral of the form 



F (&S dx can be shown to have 
\dx/ 


a stationary value for ^ = Constant (reference k, for example). Thus, 


dx 

since the integrand of equation (l) is a function only of the slope 
the profile of minimum drag has plane surfaces. Then, with -a = ^ and 
the notation of figure 1, the drag coefficient may be written as 


c d = ^ iVi + 2 ( 1 - x i ) iS 2 e> 2 - yb ^ 


( 2 ) 


where x-^ is the location of maximum thickness and 


s P K ) 




Then with 


% - s : 


0 t h 

- -.2 rrw 


h = 1 - 


^b 


where t is the thickness coefficient, equation (2) may be written as 


= t P-l - hP 2 + (h - lJP-hj 


(3) 
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Reference 3 shows that 



w 


where M ro is the stream Mach number and 7 is the ratio of specific 
heats at constant pressure and constant volume, is in excellent agree- 
ment with the pressure coefficients given by the exact shock or expan- 
sion relations for 1.5 - Mo, < 00 • It may be noted that this expression 
approaches both the linearized value and the hypersonic value in the appro- 
priate Mach number ranges. With 0 = md as the independent variable, 
equation (!<-) becomes 

P = K(0) 

m^c 

where 



Then, with these relations substituted into equation ( 3 ), the drag 
coefficient becomes 


= f-L - hK 2 + (h - 1)| m^pTl 

m^c L ^ _j 


(5) 


where 


= 


mt 


mt h 


1 2x n 


0 2 2 1 - x. 


K l sK ( e l) K 2 = k( 0 2 ) 
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Values of and h are required which make the drag a minimum. 

Thus, solutions are required for the equations 


Sc c 

dx n 


= 0 


dc d 

dh 


= 0 


for 


*1 


and h. 


The hase pressure coefficient is dependent upon Mach number, airfoil 
shape, Reynolds number, the nature of the boundary layer, and so forth. 

Of these the Mach number is by far the most' important parameter. Since 
the base pressure depends to some extent on the base height, a term 


involving 


* P b 

dh 


should be included in the second equation used to deter- 


mine the airfoil of minimum drag. The omission of this term however is 
of little consequence. Its effect is to modify slightly the effective 
value of the base pressure. For turbulent boundary layers at high 
Reynolds numbers this change is small and the change in airfoil shape 
due to the neglect of this term is of little consequence. If the term 


Sp b 

dh 


is omitted, the equations defining the airfoil of minimum drag are 


dXn 


k 

m^c 


e 




+ e 2 %' 


= 0 


( 6 ) 


dc d 

dh 


2t_ 

m p c 


(*2 


e 2 Kg’ 


+ 1 * 2 p-b 


= 0 


(7) 


where 


K 1 - 2K = £? 

30 e^/e 2 + c 2 


From the solutions of equations (6) and ( 7 ) for 0-^ and 9 2 as 
functions of (m^ 2 - l)^, the location of maximum thickness and the 
trailing-edge thickness are determined from 
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*1 = 


mt 

2 ©n 




( 8 ) 


t 


= 1 - h = 


1 + 0 2 mt 


_1_ 


0 1 0 2 < < \ 
■ 2 ” t y 


The lower limit for mt insures that the rear upper and lower surfaces 
can intersect only at the trailing edge. The upper limit for mt 
insures that the location of maximum thickness will lie between the 
leading and trailing edges of the airfoil. 

It is important to note that 0-^ and 0^ are functions only of 
(Moo^ - l)Pb* Therefore, and y^/t are functions only of t - 1 

and (mJ^ - ljPb* The number of parameters has thus been reduced from 
three to two. From equation (8) the trailing-edge thickness is zero 

0-,0p i 

for mt = -2- — ; thus, the profile of minimum drag has a sharp 

01 - 0 2 

trailing edge for this condition. For profiles with a sharp trailing 
edge, only equation (6) need be solved and the location of maximum 

thickness then is a function only of tVM^f - 1. For mt = 20^, the 
the position of maximum thickness for minimum drag is at the trailing 
edge. For values of mt > 20 no mathematical minimum for the drag 
exists. 

The solution of the equations defining the profile of minimum 
pressure drag is especially simple for linearized flow (reference 2, 
for example) . The approximate expression for the pressure coefficient 
given by linearized theory- is 


P = 


|4- 2 


M_ - 1 


and the equations which determine the profile of minimum drag are 


I * - M '* 2 - 2 + = 0 or 9l = ~ 92 


(9) 
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^ + A>) ' 0 or 9 2 = X (10) 

Then, the location of maximum thickness and the trailing-edge thickness 
axe found to he 


2t 

X 1 “ mP-K 




mP, 


t “ 2t + 2 


2 - — 
x 


A, < . < A 

-r- s-/ 


For a profile with a sharp trailing edge 


(mt < - -A), 


(H) 


the 


familiar re suit that, for minim um drag, the location of maximum thickness 
is at the midchord is found as the solution of equation (9) for y b = 0. 


The solution of equations (6) and (7) is shown in figure 2 where 0p 
and 0 2 are given as functions of (m * 2 - l)P b . The dashed vertical 

line at -(m ^ 2 - l)Pb equal to 1.428 is the limiting value of the 
vacuum pressure coefficient for infinite Mach number. The solutions for 
linearized flow are also shown in figure 2 for comparison. It may he 
noted that the values of 0p for linearized flow are in good agreement 
with the results of this analysis, whereas 0 2 shows appreciable devia- 
tion. Curves of t / m ^ 2 - 1 are also shown as functions of 

in figure 2. The upper set of curves gives the value of t^M^ - 1 for 
which the opt imum profile has the maximum-thickness location at the 

trailing edge and the lower set of curves gives the value of ty - 1 
for which the optimum profile has a sharp trailing edge. At a given value 

of (mJ^ - l)p b , true minimum-drag profiles do not exist for values of 

t^Moo 2 - 1 greater than the values given hy the upper curve. The profile 
of least drag, however, under those conditions is one which has the 

maximum -thickness location at the trailing edge. For values of t^Mco 2 - 1 
less than those of the lower curve, the profile of minimum drag has a 
sharp trail ing edge an d the position of maximum thickness for minimum drag 
may he determined from figure 3. Profiles for minimum pressure drag are 
readily determined with the use of figures 2 and 3- For a given base 
pressure coefficient and Mach number, values of 0 b and 0 2 are found 


(Mo 2 - ljPb 
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from figure 2. 


Then for 


“irO < / P < 

- 2g -q = tyMa, - 1 = 2 & 2 _> the position 


of maximum thickness and the trailing-edge thickness are found from 
equations (8) for a given thickness coefficient . For values of 



^ 1^2 

■2 3 2 —, the profile has a sharp trailing edge and the 

e l ' d 2 


position of maximum thickness for minimum drag is determined directly 
from figure 3 • 


In order to show more clearly the variation of the shape of the 
airfoil for minimum drag with Mach number, a number of profiles have 
been calculated for the base pressures shown in figure 4. The base 
pressure coefficients in this figure are, in part, based upon some 
knowledge of the base pressures for turbulent boundary layers and, in 
part, upon the variation of the vacuum pressure coefficient with Mach 
number. For comparison the vacuum pressure coefficient By is also 
shown. 


In figure 5 the position of maximum thickness and the trailing- 
pdge thickness for minimum drag are presented as functions of Mach 
number for the base pressure coefficients of figure 4. Here it may be 
noted that, for a given thickness ratio, the minimum-drag airfoil has 
a sharp trailing edge for the lower Mach numbers and at higher Mach 
numbers the trailing edge is blunt. Further, the Mach number at which 
the minimum-drag profile first has a blunt trailing edge is higher for 
the thinner airfoils. For a given thickness ratio, the position of 
maximum thickness moves rearward with increas in g Mach number until it 
is located at the trailing edge. 

Also shown in figure 5 are the optimum profiles for the 6 - and 
10 -percent -thick airfoils calculated from the linearized equations (ll). 
For the 6 -percent-thick airfoil the position of maximum thickness remains 
fixed at the midchord up to a stream Mach number of approximately 5* At 
Mach numbers greater than the position of maximum thickness moves 
rearward and the trailing-edge thickness increases in a manner similar 
to that determined from the nonlinear equations. For the 10-percent- 
thick profile the position of maximum thickness and the trailing-edge 
thickness show an erratic behavior at the low Mach numbers, that is, 
at low values of the base pressure coefficient. At a stream Mach number 
of 1-5 the optimum position of maximum thickness is at the O .56 chord. 

At Mach numbers from 1.5 to 2.1 the maximum-thickness location moves 
forward, and at higher Mach numbers it moves rearward in a manner similar 
to that determined from the nonlinear equations. In general, the linear- 
ized theory predicts the position of maximum thickness fairly accurately 
over a wide Mach number range but does not predict the optimum trailing- 
edge thickness so well. 
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In order to illustrate the effect of the airfoil geometry on drag, 
the drag coefficients of the optimum profiles have teen calculated for 
the 6- and 10-percent-thick airfoils. In figure 6 the ratio of the 
optimum drag coefficient c. to the drag coefficient of a double- 

a opt 

wedge airfoil c^ is shown as a function of the stream Mach number 

for the profile determ i ned by both the nonlinear and linear theory. For 
the 10-percent-thick airfoil the drag presented is not a TnlnlTnrrm above 
a Mach number of approximately 6. The drag curve above this Mach number 
Is that for an airfoil with the maximum thickness located at the trailing 
edge. Here it may be noted that the optimum profiles have substantially 
lower drag than the double-wedge profiles, particularly at high Mach 
numbers. Further,, the drag reduction is considerably greater for the 

10- percent-thick airfoil than for the 6 -percent -thick airfoil. The 
dashed' curves of figure 6 correspond to the drag coefficients for the 
profiles found by linear theory. In calculating the drag, however, the 
nonlinear expression for the pressure coefficient was used; that is, 
only the geometry of the airfoils was calculated by linear theory. Here 
it may be noted that the drag for the profiles determined by linear 
equations is only slightly higher than the drag of the optimum profiles 
found with the nonlinear relations. The largest difference in drag 
between the 6-percent -thick profiles determined from the nonlinear equa- 
tions and the profiles determined from the linear equations is less than 

11- percent; for the 10-percent-thick profiles the differences in drag are 
even less. 

Profile for minimum drag for a given area . - The problem of determining 
the profile of minimum drag and satisfying a given structural condition 
is an isoperimetric problem of the calculus of variations. The equations 
for determ ini ng the supersonic profile of minimum drag for a given struc- 
tural, condition have been developed in reference 2. The equations for 
determining the airfoil of minimum drag for a given area (or torsional 
strength) are 





y(0) = 0 
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where y' = X is a constant to he determined from the solution, 

A is the area, and the subscript HE refers to the values at the 
trailing edge of the airfoil. The solution of these equations for the 
linearized form of the pressure coefficient gives the profile shape as 



where x-^ is the position of maximum thickness and is given by 

mP, 

tzm 

' vr 

and the thickness coefficient is given by 



The solution of equations (12) for the nonlinear form of the pressure 
coefficient was obtained by an iterative procedure. 

Figure 7 presents a comparison of the profiles determined by linear 
and nonlinear theory for A = 0.05 for several Mach numbers and for 
the base pressure coefficients given in figure 4. The linearized theory 
gives a location of maximum thickness farther forward and a smaller 
tralling-edge thickness and thickness ratio than given by nonlinear theory. 
These differences in geometry of the profiles determined from linear 
and nonlinear theory follow the same trend as for the airfoils' of a given 
thickness ratio; this trend may be expected for other auxiliary structural 
conditions as well. The drag of the profiles found by linear and non- 
linear theory (equation (4) was used for computing the pressures for the 
evaluation of the drag in each case) differed by less than 2 percent for 
Mach numbers from 2 to 10. 


CONCLUDING REMARKS 


The results show that the profile shape for minimum drag for a given 
thickness ratio or foh a given area is determined with sufficient accuracy 
by linear theory over the entire supersonic Mach number range since the '■drag 
coefficients for these profiles are only slightly higher than those for 
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optimum profiles determined "by nonlinear theory. It would appear that 
linear theory should also he adequate for dete rminin g profiles for 
minimum drag for other auxiliary structural conditions since moderate 
deviations from the opt imum shape have only a sm al l influence on the 
pressure drag coefficient. 

It appears that, when airfoils with finite trailing-edge thicknesses 
are considered, the linearized theory may he used for deter minin g profiles 
of mi nirmim drag (although not the drag itself) at least up to Mach numbers 
of 10. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va., November l4, 1951 
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Figure 1.- Optimum profile for a given thlcimesB ratio. 
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Figure 2.- The parameterB 0-^ and &2 for profiles of minimum, drag as 

functions of "base pressure and the limits on t\^l ^ - 1 for which 
the solutions are valid. 
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Figure k.~ Vacuum pressure coefficient and base pressure coefficient as 

functions of stream Mach number. 
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Position of maximum thioknesB, Traillng-edge thickness 
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Figure 5 .- Position of maximum thickness and trailing-edge thickness for 
airfoils of minimum drag for base pressures presented in figure b. 
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Figure 7»- Cong>arison of optimum profiles for a given area as determined 
"by linear and nonlinear theory. A = 0.05. 
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